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Abstract 

The  stress  and  strain  distribution  in  an  impulsively  blastloaded  plate  is  far  from  equili¬ 
brium.  Part  of  the  plate  will  be  elastically  strained  whilst  other  part  will  be  strongly  pla¬ 
stically  deformed.  A  high  explosive  (HE)  charge,  detonating  near  to  a  flat  steelplate  pro¬ 
duces  an  impulsive  aerea-load.  Stress  and  strain  is  time  dependend  and  is  different  from 
element  to  element  in  the  plate.  The  time-dependend  displacement  of  nodes  in  the  net 
of  finite  element  calculation  and  the  time  dependend  stress  and  strain  in  elements  that 
correspond  to  the  nodes  will  be  discussed. 
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1.  Introduction 

The  designer  of  blast  resistant  structures  should  be  aware  that  in  an  impulsively  blastloa¬ 
ded  plate  the  distribution  of  stress  and  strain  is  far  from  equilibrium.  Static  design  prin¬ 
ciples  are  based  on  equilibrium  and  must  be  used  carefully. 

A  high  explosive  (HE)  charge,  detonating  near  to  a  plate  produces  an  impulsive 
aerea-load.  An  impulse  is  imparted  into  the  plate  and  any  element  of  the  plate  is  set 
into  motion.  Only  a  small  fraction,  in  the  order  of  1%  of  the  kinetic  energy,  can  be 
transformed  into  elastic  deformation.  Most  of  the  kinetic  energy  must  be  transformed 
into  plastic  deformation.  Stress  and  strain  is  time  dependend  and  is  different  from  ele¬ 
ment  to  element  in  the  plate.  An  extreme  non-equilibrium  loading  situation  was  chosen 
in  order  to  highlight  the  differences  to  static  loading. 

The  time-dependend  displacement  of  nodes  in  the  net  of  finite  element  calculation  and 
the  time  dependend  stress  and  strain  in  elements  that  correspond  to  the  nodes  will  be 
discussed. 
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2.  Test  Arrangement,  Experiments 

A  spherical  HE-charge  (1  kg  of  PFTN)  detonates  at  some  distance  (HOB)  above  the 
center  of  a  square  steelplate  (1*1*0.002  m3).  The  plate  was  clamped  along  all  edges  to 
a  rigid  support.  The  test  arrangement  is  shown  in  Figs. 2.1  and  2.2  (Ref.  1,  2). 

A  typical  experimental  result  is  shown  in  Figs.2.3  and  2.4.  The  final  shape  of  the 
impulsively  loaded  plate  is  close  to  a  pyramid,  whilst  a  statically  loaded  plate  results  in 
a  near  spherical  shape.  If  the  final  shape  is  different,  there  must  be  a  different  mecha¬ 
nism  of  deformation. 

The  maximum  deflection  of  the  plate's  centerpoint  as  a  function  of  the  specific  blastim- 
puls  (impulse  per  unit  area  imparted  to  the  plate)  is  shown  in  Fig.2.5.  The  results  from 
10  experiments  with  a  plate  lm  *  lm  and  2  mm  thick  and  HE-charges  at  different  HOB 
and  8  results  from  numerical  calculations  obviously  lie  at  one  straight  line.  This  dia¬ 
gram  demonstrates  that  the  DYNA  3D  FE-code  is  suited  to  describe  the  large  deformati¬ 
on  of  a  steelplate  under  impulsive  loading.  More  checks  were  done  to  confirm  this  re¬ 
sult. 


3.  Numerical  Calculation 

The  overall  deformation  of  the  steelplate,  as  well  as  stress  and  strain  in  different  ele¬ 
ments  of  the  steelplate  was  calculated  by  means  of  the  DYNA  3D  code.  Details,  how  to 
find  the  best  net-discretisation,  to  minimize  the  machine  time  and  to  optimize  the  input- 
and  output  procedures  can  be  found  in  the  references  3,  4,  5,  6. 

The  numerical  work  with  the  FE-code  can  only  be  successful,  if  two  sets  of  input  para¬ 
meters  are  correct: 


-  the  load  function  in  spacial  and  temporal  distribution. 

-  the  material  properties. 

The  input  load  function  was  measured  carefully  with  piezoelectric  pressure  gages  at  a 
nonresponding  platform  (Ref.l).  At  close-in  detonation  the  spherical  shockfront  is  not 
flat  when  it  impinges  the  flat  plate.  It  impinges  first  at  the  center  of  the  plate  (ground 
zero)  and  then  spreads  to  the  edges.  Different  pressure- time  histories  for  up  to  20  zones 
at  the  plate's  surface  were  determined.  I 

Additional  numerical  calculation  was  done  for  arrangements  where  no  direct  compari¬ 
son  with  experimental  results  was  planned.  Identical  input  blast  parameters  had  to  be 
used  for  different  arrangements.  In  that  case  blast  parameters  from  the  literarature  were 
taken  (Ref.5;  6). 


Most  of  the  numerical  calculation  was  done  with  standard  values  for  the  material  proper¬ 
ties  of  mild  steel  under  "static"  load: 


Mass  Density 
Modulus  of  Elasticity 
Poisson's  Ratio 
Yield  Stress  Tension 
Yield  Stress  Shear 
Tangent  Modulus 
Temperature 


7770  kg  /  m3 

207  GPa 

0.3 

340  MPa 

195  MPa  (0.6  of  tension) 

68.9  MPa  /  Strain  (m/m) 

20  degree  C 
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Numerical  calculation  was  done  in  the  Held  of: 

Different  HoB  for  1  kg  HE  charges  (Ref.l,  4) 

Different  spatial  and  temporal  pressure  distribution  at  the  flat  plate  for  spherical  blast 
(Ref.  3,  4) 

Centric  and  excentric  detonation  (Ref.  3) 

Square  and  rectangular  plates  (Ref.l;  3) 

Effect  of  tension  yield  stress  (525  MPA  instead  of  340  MPa)  (Ref.  1) 

Effect  of  tangent  modulus  (3  GPa  instead  of  69.9  MPa)  (Ref.  1) 

Effect  of  stiffeners  (Ref.  4) 

Effect  of  a  girder  (Ref.5) 

Discretisation  of  the  plate  in  2  layers 

Results  from  different  reports  will  be  discussed. 


4.  Phenomenology  of  Impulsive  Blastloading 

A  high  explosive  (HE)  charge,  detonating  near  a  plate  produces  an  impulsive  aerea- 
load.  The  time-scale  of  the  process  generally  is  milliseconds.  The  reflected  blastpressu- 
re  at  the  plates  surface  rises  suddenly  to  its  peak  overpressure  and  decreases  to  zero  in  a 
time  that  is  short  relatively  to  the  reaction-time  of  the  plate. 

An  impulse  is  imparted  into  the  plate  and  any  element  of  the  plate  is  set  into  motion. 
The  plate's  mass  contains  an  amount  of  kinetic  energy  that  must  be  transformed  into  de¬ 
formation.  Only  a  small  fraction,  in  the  order  of  1  %  of  the  kinetic  energy,  can  be  trans¬ 
formed  into  elastic  deformation.  Most  of  the  imparted  energy  must  be  transformed  into 
plastic  deformation.  The  reaction-time  is  over  when  all  the  kinetic  energy  is  transfor¬ 
med  and  the  plate  is  at  rest. 

The  acceleration,  velocity  and  displacement  vs.  time  of  a  selected  point  (node)  in  the 
center  of  the  square  plate  will  be  discussed.  The  shape  of  the  plate  at  different  moments 
of  time  will  be  shown  as  well  as  the  final  shape  after  the  loading. 

Fig.4.1  The  shape  of  the  plates  cross-section  from  the  center  (left)  to  the  fixed  edge  in 
the  middle  of  the  span  at  500  mm  (right)  is  shown  at  4  moments  of  time  after  the  shock- 
front  arrival  (Ref.2).  The  plate  was  loaded  from  above  and  moves  downwards.  At 
0.125  ms  the  movement  starts  with  a  flat  plate.  At  0.875  ms  the  flat  bottom  has  rea¬ 
ched  a  displacement  of  30  mm.  A  "knee"  moves  along  the  plate  from  the  fixed  edge  to 
the  center  and  causes  a  rim  at  constant  slope.  The  section  of  the  plate  that  was  passed 
by  the  knee  is  at  rest.  The  flat  bottom  moves  downwards  at  constant  velocity.  It  gets 
smaller  with  time.  At  2.2  ms  the  knee  has  reached  the  plate's  center.  The  whole  plate 
was  plastically  deformed  and  is  at  rest  now.  This  is  the  end  of  the  reaction  time. 

Fig.4.2  The  final  cross-sectional  shape  of  the  deformed  steelplate  for  different  blastim- 
pulses  is  shown.  The  plate  was  fixed  at  the  left  side.  Different  impulsive  load  was  pro¬ 
duced  by  different  distances  (HoB)  of  the  HE-charge  above  the  plate.  The  slope  of  the 
rim  gets  steeper  with  increasing  impulse.  The  final  shape  was  reached  at  the  same  time 
of  2.2  ms,  independently  of  the  load. 

Fig.4.3  The  center  of  the  square  plate  was  accelerated  during  the  positive  pressure 
duration  of  0.22  ms.  Maximum  acceleration  is  58,000  g's.  Acceleration  is  zero  from 
0.22  ms  to  1.8  ms,  when  a  force  acts  in  the  opposite  direction.  Maximum  acceleration 
in  the  opposite  direction  is  36,000  g's  at  2.2  ms. 


163 


Fig.4.4  A  certain  amount  of  kinetic  energy  was  impaled  into  the  plate  during  the  phase 
of  acceleration.  We  are  intersted  in  the  mechanism  that  reduces  thd  kinetic  energy.  At 
2.2  ms,  at  the  end  of  the  reaction  time,  there  is  no  kinetic  energy  left,  the  plate  is  at 
rest. 

Fig.4.5  The  velocity  reaches  a  constant  value  of  75  m/s  at  the  end  of  the  acceleration 
phase.  The  plate's  centerpoint  moves  at  constant  velocity  from  0.22  ms  to  1.8  ms  and 
then  is  rapidly  stopped.  Zero  velocity  is  reached  at  2.2  ms. 

Fig.4.6  The  centerpoint  moves  at  nearly  constant  velocity  for  2.2  ms  where  it  reaches 
its  maximum  displacement  of  150  mm.  Some  relatively  small  vibrations  occur  later 
than  2.2  ms. 


The  discussion  of  the  diagrams  results  in  the  following  conception  of  the  deformati¬ 
on  mechanism: 

The  deformation  of  steelplates  at  impulsive  blastloading  happens  in  the  timeframe  of 
milliseconds.  An  amount  of  kinetic  energy  is  imparted  into  the  plate  and  must  be  trans¬ 
formed  into  plastic  deformation. 

Any  point  of  the  plate  is  accelerated  during  the  time  of  positive  pressure  duration. 
Extremely  high  values  of  acceleration  occur  (some  10,000  g's).  At  the  end  of  the  acce¬ 
leration  phase  a  constant  velocity  is  reached  (some  10  m/s).  A  plastic  wave  or  knee 
starts  in  the  moment  of  loading  from  all  fixed  edges  and  runs  at  constant  velocity  along 
the  plate  to  the  plate's  center  (225  m/s).  Any  point  that  is  caught  by  the  knee  will  be  re¬ 
tarded  roughly  and  comes  to  rest.  Most  of  the  kinetic  energy  must  be  transformed  in 
the  running  plastic  knee.  The  reaction-time  is  determined  by  the  time  that  it  takes  for 
the  plastic  wave  to  run  through  the  material  from  any  fixed  edge  to  the  center  of  the 
plate. 

The  mechanism  results  in  a  rim  at  constant  slope,  that  depends  on  the  loading  impulse, 
the  mass  and  the  strength  of  the  material.  The  final  shape  of  the  deformed  plate  is  pyra¬ 
midal. 

The  mechanism  of  energy  transformation  and  material  deformation  results  in  a  spatial 
and  temporal  stress  and  strain  distribution  in  the  plate  that  will  be  discussed  in  the  follo¬ 
wing  section. 


5.  Stress  and  Strain  in  the  Plate 

The  effect  of  a  blastwave  on  a  steelplate  was  computed  by  means  of  the  explicit,  3-di¬ 
mensional  FE-code  DYNA  3D  (Ref.5,  6).  Acceleration,  velocity  and  displacement  at 
different  points  (nodes)  of  the  plate  was  calculated  as  well  as  the  shear  stress,  the  effec¬ 
tive  stress  and  the  effective  plastic  strain  in  corresponding  elements. 

A  spherical  150  kg  TNT-charge  detonates  at  a  distance  of  6  m  from  a  steelplate.  The 
plates  dimensions  are  7.2m  *  2.5m  *  0.004m.  The  blast  data  are  3.66  MPa  reflected 
peak  overpressure,  3  ms  positive  duration  and  2564  Pa*s  specific  reflected  blastimpul- 
se.  The  yield  stress  was  340  MPa.  Large  plastic  deformation  occurs. 

The  time-scale  is  10  ms  in  all  diagrams  shown  in  this  section. 
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Fig.5.1  A  pressure  pulse  of  3660  kPa  peak  overpressure  and  3  ms  positive  duration  im¬ 
pinges  on  the  flat  surface  of  the  plate  at  time  0.  It  transfers  impulse  to  the  plate  and  in¬ 
duces  the  plate  to  movement.  According  to  its  mass  and  velocity  the  plate  contains 
kinetic  energy  that  has  its  maximum  at  1.8  ms.  The  kinetic  energy  decreases  continous- 
ly  and  is  near  0  at  6  ms. 

Fig.5.2  The  effective  stress  in  3  elements  of  the  plate  is  shown  as  a  function  of  time. 
The  element  B  is  situated  near  to  the  clamped  edge  and  the  element  J  at  the  midspan. 

It  takes  0.2  ms  to  reach  the  yield  stress  near  the  edge  (B)  and  1ms  at  the  midspan  (J). 
The  stress  rises  very  fast  (sub-millisecond)  in  each  element. 

The  sudden  rise  of  stress  starts  at  different  time  in  the  elements  from  "near  zero"  stress 
to  yield  stress.  The  "mechanism"  that  causes  the  stress  to  reach  yield  stress  runs  with 
longitudinal  acoustic  wave  velocity  (5800  m/s)  from  the  clamped  edge  to  the  plate's 
center. 

In  all  three  elements  the  stress  remains  constant  at  a  "yield  stress  plateau"  for  a  period 
of  time.  The  moment  of  time  is  marked  by  an  arrow  when  the  stress  in  an  element 
exceeds  the  yield  stress.  This  happens  first  at  the  clamped  edge  (B)  and  2ms  later  at  the 
midspan  (J).  There  must  be  a  "mechanism"  that  causes  the  stress  to  exceed  yield 
stress  at  different  moment  of  time  in  different  elements. 

When  the  Kinetic  energy  returns  to  zero  (KE  =  0  at  6  ms)  the  effective  stress  in  the  ele¬ 
ments  decreases  below  yield  stress. 

Fig.5.3  The  time  history  of  the  shear  stress  in  the  elements  is  nearly  identical  to  the  ef¬ 
fective  stress.  Shear  yield  stress  was  taken  as  0.6  times  the  tension  yield  stress.  Shear 
stress  is  an  important  parameter  in  the  deformation  mechanism  of  the  steelplate. 

Fig.5.4  Effective  plastic  strain  arises  in  an  element  at  the  time  when  the  stress  exceeds 
yield  stress  in  this  element.  The  time  to  reach  the  maximum  strain  is  different  in  the 
elements. 

The  strain  rate  in  m/m/s  can  be  calculated  and  is  different  in  the  elements. 

When  the  Kinetic  energy  returns  to  zero  (KE  =  0  at  6  ms)  the  plastic  strain  in  the 
elements  has  reached  its  maximum  and  constant  value. 

Fig.5.5  and  5.6  Stress  and  strain  in  4  elements  along  the  short  fixed  span  are  shown. 
The  strain  exceeds  the  yield  stress  at  the  same  moment  in  all  the  elements.  Highest 
stress  arises  in  the  middle  of  the  span  (A).  In  the  comer  (I)  the  stress  exceeds  yield 
stress  marginally.  The  plastic  strain  is  marginal  in  the  comer  (I,  0.25%)  and  maximum 
in  the  middle  (A,  4.2%).  Strain  at  the  short  span  is  smaller  than  at  the  long  span  (see 
Fig.5.4,  Element  B). 

Fig.5.7  The  maximum  strain  in  an  element  is  directly  proportional  to  the  maximum 
stress  in  the  element.  The  stress-time  history  does  not  have  an  effect  on  the  final  plastic 
strain. 
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Fig.  1  Test  Arrangement: 

Clamped  Quadratic  Steel  Plate 
Parameters : 

Plate  Side  Length  l 
Plate  Thickness  h 
Charge  Height  HoB 


i 


Fig.  2  The  lest  arrangement  -  a  flat  square 

steelplate  lm*1m  clamped  and  bolted 
to  a  rigid  support.  The  1kg  charge  made  of 
plastic  explosive  is  wrapped  into  thin  fabric 
and  detonated  at  some  height  above  the 
plate  center. 


CT\ 

CO 


10  Experiment's 
5  ADINA 
3  DYNA  3  D 


(Ref  .1) 
(Ref  .1) 
(Ref. 3) 
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Fig.  3  Typical  Shape  of  Homogeneous  Fig.  4  A  clamped  square  steelplate  deforms  to  a 

Steel  Plates  at  Static  and  pyramid  at  high  impulsive  loading. 

Impulsive  Loading  The  straight  line  pattern  at  the  flat  plate 

remains  straight  at  the  pyramid  after  deformation. 


t=  0.125  ms 


t=  0.875  ms 


Fig. 4.1  The  displacement  of  the  plates  cross-section 
from  the  centerpoint  to  one  clamped  side  at 
different  instants  of  motion. 
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Fig. 4. 2  Final  Shape  of  the  Deformed  Steelplate  at  Different  Blastload 

Square  Steelplate  1x1  m2;  0.002  m  thick 
Final  Position  at  2.2  ms;  Plastic  Wave  Velocity  225  ms"^ 
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Fig. 4. 6  Displacement  vs  Time 

Centerpoint  of  a  Blastloaded  Square  Steelplate 


Effective  Stress  (MPa) 


Fig.5.1  Pressure  Pulse  and  Time  History  of  the  Kinetic  Energy 
in  the  Plate 


KE=0 


Fig.5.2  Effective  Stress  in  3  Elements 
Element  B  -  near  the  edge 
Element  J  -  at  midspan 
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Effective  Plastic  Strain  (%)  ~  Maximum  Shear  Stress  (MPa) 


Shear  Yield  Stress 
is  taken  equal  to  0.575 
the  tension  yield  stress 
in  DYNA 
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Fig. 5. 5  Effective  Stress  vs  Time 

Four  Elements  along  the  clamped  edge 


Fig. 5.6  Effective  Plastic  Strain  (%)  vs  Time 
Four  Elements  along  the  clamped  edge 
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Fig. 5. 7  stress  exceeding  yield 

stress  vs  permanent  strain  (DYNA) 


